Chiral Symmetry Breaking and Scalar Confinement 



P. BicudcE and G. M. MarquetQ 
Dep. Fisica and CFIF, Instituto Superior Tecnico, Av. Rovisco Pais 1049-001 Lishoa, Portugal 

We address the old difficulty in accommodating the scalar quark-antiquark confining potential 
together with chiral symmetry breaking. We develop a quark confining potential inspired in the 
QCD scalar flux tube. The coupling to quarks consists in a double vector vertex. We study the 
Dirac and spin structure of this potential. In the limit of massless quarks the quark vertex is vector. 
Nevertheless symmetry breaking generates a new scalar quark vertex. In the heavy quark limit the 
coupling is mostly scalar. We solve the mass gap equation and find that this potential produces 
spontaneous chiral symmetry breaking for hght quarks. The quantitative results of this model are 
encouraging. 
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I. INTRODUCTION 

Spontaneous Chiral Symmetry Breaking (SxSB) is ac- 
cepted to occur in low energy hadronic physics. Another 
important feature of hadronic physics, suggested by the 
spectroscopy of hadrons, by lattice simulations and by 
models of confinement is scalar confinement. However 
SxSB and scalar confinement are apparently conflicting, 
since the first requires a chiral invariant coupling to the 
quarks, like the vector coupling of QCD. We address a 
recent quest of Bjorken 1], "how are the many disparate 
methods of describing hadrons which are now in use re- 
lated to each other and to first principles of QCD?" . Al- 
though the vector confinement of quarks is not yet ruled 
out here we try to solve this old chiral symme- 

try versus scalar confinement conflict of hadronic physics, 
which remained open for many years. 

In this paper we explore scalar conflnement from the 
perspective of chiral symmetry breaking. In SectionlTllwe 
motivate the importance of both SxSB and scalar con- 
finement. We show these features of hadronic physics to 
have some subtle weaknesses that we capitalize to con- 
struct a model. The potential used in our Quark Model 
(QM) is defined in Section IIIII The self consistent mass 
gap equation for the quarks is derived in Section llVI In 
Section we solve numerically the mass gap equation 
and calculate the quark condensate. Finally, in Section 
IVllwe present some conclusions. 



II. MATCHING CHIRAL SYMMETRY 
BREAKING WITH SCALAR CONFINEMENT 

The QCD lagrangian is chiral invariant in the limit 
of vanishing quark masses. Nambu and Jona-Lasinio 
showed that including chiral symmetry in fermionic sys- 
tems provides a natural explanation for the small pion 
mass, which is much lighter than all the other isovector 
hadrons. Because of this crucial fact the mechanism of 
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SxSB is accepted to occur in low energy hadronic physics 
for the light flavors u, d and s, where mmmd ^ TOs < 
^QCD < Mn /i. Similarly to the vector Ward identities 
in gauge symmetry, the axial Ward identities constitute a 
powerful tool of chiral symmetry. The techniques of cur- 
rent algebra led to beautifully correct theorems, the Par- 
tially Conserved Axial Current (PCAC) theorems. The 
different variants of QMs are widely used as simplifica- 
tions of QCD. They are convenient to study quark bound 
states and hadron scattering. Recently we have shown 
these beautiful PCAC theorems, like the Weinberg theo- 
rem for TT — TT scattering, to be reproduced by QMs with 
SxSB. Another important benefit of having SxSB in the 
QM is the reduced number of parameters. The mass 
gap equation generates a dynamical constituent quark 
mass, which is no longer an independent parameter, even 
for quarks with a vanishing current mass. The quark- 
quark, quark-antiquark, antiquark-antiquark potentials, 
and the quark-antiquark annihilation and creation in- 
teractions are all originated in the same chiral invari- 
ant Bethe-Salpeter kernel. Therefore any QM for light 
quarks should comply with the SxSB. Moreover the mi- 
croscopic coupling of the quark to the confining interac- 
tion should include the vector coupling which is present 
in the quark-gluon vertex of QCD. However there is some 
evidence that a vector quark-quark potential is not suf- 
ficient to provide the expected scale of SxSB of the or- 
der of 200 — 300 MeV which is present both in the con- 
stituent quark mass and in the quark condensate. It was 
realized by Adler that a linear confinement with vec- 
tor couplings was not sufficient to provide the correct 
quark condensate. Moreover the gluon propagator ex- 
tracted from the lattice, when used in a one gluon ex- 
change truncation of the quark mass gap equation, is not 
able to provide the expected quark condensate 8]. This 
also happens with the gluon propagator extracted from 
the solution of truncated Schwinger-Dyson equations of 
QCD 1^. Importantly, these gluon propagators exhibit 
a non-perturbative mass. This mass should produce a 
Meissner effect in Yang-Mills fields, and this is expected 
to produce a confining string for the quarks. Here we will 
estimate the effect of the confining string on the quark 
condensate. 

On the other hand the confining potential for con- 
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stituent quarks is probably scalar. We can learn much 
by comparing simply the spectrum of the hydrogen atom 
with the masses of all hadron families. In a perturba- 
tive QCD scenario, the hadron spectroscopy would be 
dominated by the one gluon exchange, which is quali- 
tatively similar to the one photon exchange interaction 
that explains in detail atomic physics. It turns out that 
all the hadronic families, say of mesons or baryons, with 
light or heavy flavors, show similar differences with the 
hydrogen spectrum. It is remarkable that the Spin- 
Orbit potential (also called fine interaction in atomic 
physics) turns out to be suppressed in hadronic spec- 
tra since it is smaller than the Spin-Spin potential (hy- 
perfine interaction). This constitutes an evidence of non- 
perturbative QCD. Another evidence of non-perturbative 
QCD is present in the angular and radial excitations of 
hadrons, which fit linear trajectories in Regge plots, and 
suggest a long range, probably linear, confining potential 
for the quarks. This led Henriques, Kellett and Moor- 
house to develop a QM where a short range vector 
potential and a long range scalar potential partly cancel 
the Spin-Orbit contribution. The short range potential 
is Coulomb-like (inspired in the one gluon exchange) and 
its quark vertex has a vector coupling structure ipj^^ip. 
The long range pote ntial has scalar coupling ipip and is a 
linear potential [IS [H El El . 

The same scalar confinement picture is extracted from 
lattice simulations. In quenched lattices which simulate 
the quark-antiquark potential in the heavy quark limit, 
the pattern of spin-spin, tensor, and spin-orbit interac- 
tions is compatible with a scalar confinement 0, 0, . 
Moreover the presently favored confinement picture in 
the literature is the flux tube, or string picture, with 
tension a ~ 200 MeV/Fm. Quantum mechanics suggests 
that a thin string, in its ground-state, should be a scalar 
object E3- Only higher energy excitations of the string 
would have angular momentum. This was capitalized by 
Isgur and Paton in the flux tube model El ■ 

In this paper we will assume the quarks are coupled to 
a scalar object which provides a linear confinement. On 
top of that we want this coupling to use, as a microscopic 
building block, the vector gluon-quark which is present 
in QCD. Notice that to get a Lorentz scalar coupling 
a simple one gluon vertex is not enough. The coupling 
needs at least two vertices. The simplest way to achieve 
this is to have the string emitting two effective gluons 
that couple to the same quark line. This double vertex, 
coupling the string to a quark line via two intermediate 
gluon propagators, is presented in Fig. ^ Effectively 
this double vertex is similar to the vertices that couple 
a quark to a gluon ladder in the soft pomeron models 
|l9l I20I . It is also related to the light quark vertex in 
the heavy-light quark bound states studied in the local 
gauge coordinate 0, lU I3 ■ Such a double vertex was 
also introduced in the coupling of short strings to quarks 
|24j . In the cumulant expansion formalism int terms of 
gluon correlators pH |2^] this means the quark-antiquark 
coupling will be dominated by four gluon correlators. 
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FIG. 1: The coupling of a quark to a string with a double 
gluon vertex 



III. THE DOUBLE VERTEX 
NON-PERTURBATIVE CONFINING 
INTERACTION 

In this section we construct the simplest possible cou- 
pling, that simulates that of a quark line with a scalar 
string using two vector couplings as building blocks. The 
most general coupling of this kind is presented in Fig. Q] 
and it is simply read, 

, , (1) 

where F is the dirac structure of the fermion-gluon inter- 
action and the usual color interaction A°/2. We de- 
note the quark propagator by S{k) and the gluon propa- 
gator by G°'^{k). Finally VsW'^'^'^ represents the coupling 
of the gluon pair to the string. 

To get the coupling of the string to a light quark, we 
follow the coupling obtained in the heavy-light quark sys- 
tem, computed in the local coordinate gauge j22i]- This 
model interpolates between the heavy-light meson in the 
local gauge and the effective QM. So, for the Dirac struc- 
ture of the quark-gluon sub- vertices F we have 7*^ matri- 
ces, which is also compatible with the Coulomb gauge. 

In the color sector, as already stated, the coupling of 
the same sub- vertices has a A"/2 structure, where A" are 
the Gell-Mann matrices. The remaining sub- vertex in- 
cludes the couphng of two color octets, see Fig. El The 
string is also a colored object; it contains a flux of color- 
electric field. For a scalar coupling, which is symmetric, 
we use the symmetric structure function defined by 

{A'',A''} = -<5"^ + 2d°''^A'= . (2) 
3 

This will result in a color contribution for the effective 
vertex of 

In QMs the string usually couples to the quark line with 
a A°/2. In our case it couples with two A°/2, one for 
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8 + 8 = 27 +10 + 10 + 8 + 8 + 1 

FIG. 2: String color 

each sub- vertex. But as we can see from Q the effective 
resuh is the same. 

The gluon propagators and the different sub-vertices 
resuh in a distribution in the loop momentum k. Here 
different choices would be possible. For simplicity we 
assume that the relative momentum p — q flows equally 
in the two effective gluon lines. We also remark that 
the distribution in k is normalized to unity once the cor- 
rect string tension is included in the relative potential 
V{p — q). This amounts to consider that the momentum 
k distribution is a Dirac delta 



{27ry6^ k - 
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(4) 



If we summarize the above assumptions in a formula 
the loop vertex that will be used is 
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(5) 

The equal-time approximation, which is standard in 
QMs, allows the computation of the double vertex as a 
functional of the running quark mass nik- The quark 
mass will be computed self-consistently in the next sec- 
tions. 

We decompose the fermion propagator in the usual 
particle and antiparticle propagators 
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(6) 



ko - Ek+ ie -fco - Ek + ie 
where the quark energy projectors are 

1 + Sfe/3 + Cfe fc • a 



A+(fc) 
A-(fc) 



1 - Skf3 - Ckk a. 
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S 

Y.Vs{k)vl{k) 



(7) 



and where Sk = sin ipk — rrik / vk^ + mf^, Ck = cos cpk — 
k/ y/k'^ -\- ruk'^ and (pk is the chiral angle, a convenient 
function for algebraic and numerical computations. 
The energy loop integral can be easily calculated. 
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(8) 



J ^5(fc) = (A+-A-)/3 

= (sfc/3 + Ckk- a)P 



(9) 



Finally, using the Dirac delta distribution for the re- 
maining integrals over the three momentum k, and sum- 
ming in color indices, we get the following effective ver- 
tex, 



Veff = ^ Y (Sfc - Cfe fc • 7) 



(10) 
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In the remainder of the paper we will always assume k = 

(p + g)/2. 

Eq. H1U|) shows that the double vertex actually solves 
the problem of matching chiral symmetry breaking and 
scalar confinement. In the chiral limit of a vanishing 
quark mass, the effective vertex Voff — > —CX'^/2k ■ 7 is 
proportional to the 7'^ and is therefore chiral invariant 
as it should be, whereas in the heavy quark limit, Vcff — > 
is simply a scalar vertex. The Gell-Mann matrix 
A'^ provides the usual color vector coupling as expected 
in a QM. We anticipate that the dynamical generation of 
a quark mass will also generate a scalar coupling for light 
quarks, and this results in an effective vertex, including 
a chiral invariant vertex ^ and the standard scalar vertex 
1. 

The dependence in the relative momentum must com- 
ply with the linear confinement which is derived from the 
string picture, 



(11) 



where a ~ 200 MeV/Fm is the string constant and C 
is the algebraic color factor defined in Eq. The 
damping factor e regularizes the Fourier transform. 
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and in the limit £ — + we have 



-iVo{p-q) 



-Sttct 



(12) 
(13) 



IV. MASS GAP EQUATION 

We solve the mass gap equation using the Schwinger- 
Dyson formalism. 



s- 



(14) 



where the dressed propagator is defined in eq. (jSJ and the 
free propagator has a similar definition with the quark 
bare mass mg. We want to determine the constituent 
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FIG. 3: The self-energy term of the mass gap equation 



quark mass m,p solving the self consistent mass gap equa- 
tion. To compute the self energy in this model, using the 
effective string potential which couples to the quark line 
with double vertices, we have to deal with the three loop 
diagram of Fig. O Keeping in mind that the two dou- 
ble vertex loops are already simplified. Technically other 
diagrams, with crossed gluon legs, could exist but we as- 
sume that this is the dominant diagram. 

The mass gap equation can be much simplified in the 
spin formalism. Some useful relations we will use for this 
purpose are 

uHp) Vs'ip) = [(Tiia2)],,, 
ul{p)j3vs'{p) = Cpp ■ [(T{ia2)]ss' 
ul{p)a^vAp) - -{5'' - (1 - SpWfP) [a^{ia2)\^^, 
ul{p)(3a'vAp) - ~{sp5'' + (1 - SpWp') H(*'^2)],,, 

(15) 

With eqs . and H15() we arrive at the expected 
relation ImI 



(16) 



implying that the propagator is diagonal in the particle- 
antiparticle projection. The mass gap equation becomes 



where we used the properties of the /3 and a* matrices. 
The result for the mass gap self-energy term is 

Us{p)T.{p)vs'{p) = 

/(Pq 
J^i-^We i\p - q\) [{4 - cl)sq{^cpp)~ 



(20) 



(27r)3 

- 2skCkSq{-k + (1 - Sp)k ■ p p)+ 
+ Cqi-q + (1 - Sp)q ■ p p) + 
+ 2skCkCqk ■ q{-Cpp)- 

-2clcqk ■ q{-k + (1 - Sp)k ■ p p) ■ [cr{ia2)] 



As we can see from H18|l and (I2UII both terms of the mass 
gap equation are proportional to p ■ cr{ia2)- Since the 
Pauli matrices a are linearly independent, we can sub- 
stitute (J {102) by p and still have a mass gap condition. 
With this simplification the mass gap equation reduces 
to. 



ipSp 



(cfc - sl)sqC-p + '2skCkSqSpk ■ p- 
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- CqSpQ ■ p - 2SkCkCqCpk ■ Q- 



(21) 



2clcqSpk ■ q k p 



iV,{\p-q\)=Q. 



Notice that if we take the integrand and set q = p we 
will get X T4(0). In section we will deal numerically 
with this IR behavior. 



NUMERICAL SOLUTION OF THE MASS 
GAP EQUATION 



Usip)So ^(p)vs'(p) - Us{p)T.{p)vs'{p) = 0. (17) 

The free propagator term for a zero bare mass, mg = 0, 
is simply 

wI(p)/3(-i)^ Vs'(p) = iP SpP- [(T{ia2)],,, (18) 

The self-energy term is directly derived from the dia- 
gram of Fig. 121 using the double vertex loop VeS already 
obtained in (|10|l . 
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16 



(Sfc - Cf.)Sql3 - 2SkCkSq fe • OL + 
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(2^ 

' (Sfc + Cfe)Cg q OL + "^SkCkCq k- q(3- 

-2clcqk-qk-a) p{~i)Ve{\p-q\) , 



The mass gap equation is a difficult non-linear integral 
equation, that does not converge with the usual iterative 
methods. We developed a method to solve the mass gap 
equation with a differential equation, using a convergence 
parameter A. This parameter is the radius of a sphere 
centered init = p— qf = Oit allows us to separate the 
integral into an integral inside the sphere and another 
outside of it. 



dh 



S3 f{P,u)Veiu) . 
R3_o {2ttY 



(19) 



(22) 



In our model f{p,u) is the fmiction dependent on the 
chiral angle presented in eq. if^ . 

Let us first focus on the integral inside the sphere, 
where we have m < A. Eventually we will take the limit 
where A — > and this term will vanish. But for now we 
will expand the function / around u — and take only 
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the first non- vanishing term, 



f{p,U,Uj)Ve{u) 



27r)3 (9^2 

A(-8^g) 
87r2 



1 ^2 

duj 

-1 



Ve{u) 



ti=0 

(23) 



where uj is the cosine of the angle between p and u. For 
our particular model we have 



_ sin(2<^p) - 2p cos{2ipp)ip'p + p'^ip 



(24) 



In what concerns the integral outside the sphere (u > 
A) we can take from the beginning e = since this in- 
tegral is already regulated by A. In this case we have 



fip,u)Vo{u) 



-Sttct 
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du 



(25) 



Placing the two terms in the mass gap equation we 
finally get the equation that we can iterate to find the 
solution, 



3p Sp 



A (ip: 



2 1 

- cos{2(pp)Lp' + — sin(2<y9p)) 
p ^ p^ 
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du 



= 



(26) 



Our technique consists in starting with a large infrared 
cutoff A, where the integral term of Eq. H26|l is negligible. 
In this case Eq. H2ti|) for the chiral angle Lpp becomes es- 
sentially a differential equation which can be solved with 
the standard shooting method 27]. It turns out that 
this equation possesses several solutions, and we special- 
ize in the larger one, with no nodes, that corresponds 
to the stable vacuum '2^. Then one decreases step by 
step the A parameter, using as an initial guess for the 
evaluation of the integral the Lpp determined for the pre- 
vious value of A. In this way the integral is a simple 
function of the momentum p and we again have to solve 
a non-homogeneous differential equation. Eventually we 
are able to solve the mass gap equation for a A parameter 
which is much smaller than the scale of the interaction. 
Finally we extrapolate the set of obtained Lpp to the limit 
of A ^ 0. We test the convergence of the method com- 
puting the quark condensate , 



(27) 



The evolution of the solution as a function of the in- 
frared parameter A is clearer when we display the quark 
condensate, see Fig. ^ 

The solution of the mass gap equation is presented in 
Fig. [3 where we compare it with the single vertex model 
solution. 



VI. RESULTS AND CONCLUSION 

In this paper we build a QM for the coupling of quark 
to a scalar string. The quark confining interaction has 
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FIG. 4: Testing the convergence of the numerical method 
with the quark condensate (i/'i/'). 



a single parameter a. This QM matches the apparently 
conflicting vector coupling of QCD with a scalar confine- 
ment. Our model can be interpreted as a double vertex 
that couples the quark to the string, or alternatively as 
a 7oiS'(fc)7o vertex. Either way this vertex decomposes in 
the sum of a scalar vertex 1 and a chiral invariant ver- 
tex weighed by simple functions of the dynamical quark 
mass. In the chiral limit the scalar vertex vanishes, while 
in the heavy quark limit the confining potential is essen- 
tially scalar. Our results for the weighing factors of the 
scalar vertex and the remaining chiral invariant vertex 
are shown in Fig. 

We solve the mass gap equation for the dynamical gen- 
eration of the quark mass with the S^SB, and we in- 
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FIG. 5: The nik solutions of the mass gap equation in units 
of 




FIG. 6: Form factors for the scalar vertex (sk) and for the 
remaining chiral invariant vertex (ct) 



deed generate the constituent quark mass. We show that 
SxSB not only generates a quark mass, but generates also 
a scalar vertex for the confinement. The results are en- 
couraging because the quark condensate indeed increases 
when compared with the simpler one vertex vector con- 
fining potential. 



It is clear that the next step of this work, will consist 
in adding the shorter range one gluon exchange potential 
to the confining potential. The resulting model will have 
two parameters, one for the short range potential, and 
another one for the confining potential. These parame- 
ters will be determined in the fit of the hadron spectrum. 
In what concerns the mass gap equation, we expect that 
this will further enhance the quark condensate, possibly 
up to the expected -{230MeVf. 
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